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Abstract 

We present the full one-loop corrections to charged and CP-even neutral Higgs 
boson decays into sfermions including also the crossed channels. The calculation 
was carried out in the minimal supersymmetric extension of the Standard Model and 
we use the on-shell renormalization scheme. For the down-type sfermions, we use 
DR running fermion masses and the trilinear coupling Af as input. Furthermore, 
we present the first numerical analysis for decays according to the Supersymmetric 
Parameter Analysis project. This requires the renormalization of the whole MSSM. 
The corrections are found to be numerically stable and not negligible. 



1 Introduction 



The Higgs boson is the last not discovered particle of the Standard Model (SM) and so the 
search for the Higgs boson is the prime objective of the LHC and other future colliders. 
Apart from the SM, the Higgs boson is also predicted by its minimal supersymmetric 
extension - the Minimal Supersymmetric Standard Model (MSSM). As opposed to the 
SM, the MSSM has not only one neutral Higgs boson but it predicts the existence of two 
neutral CP-even Higgs bosons (h°, H°), one neutral CP-odd Higgs boson (A ) and two 
charged Higgs bosons (H ). The existence of a charged Higgs boson or a CP-odd neutral 
one would be clear evidence for physics beyond the SM. 

A further difference in the MSSM is the possibility for the Higgs bosons to decay not only 
into SM particles. In case the supersymmetric (SUSY) partners are not too heavy, the 
Higgs bosons can decay into SUSY particles as well (neutralinos xf, charginos xt an ^ 
sfermions f m ). The new decay channels might substantially influence the branching ratios 
of the MSSM Higgs bosons. 

At tree-level the decays into SUSY particles were studied in JTJ [2] and one-loop effects of 
the decays into charginos and neutralinos were analyzed in [3j H] and were found not to 
be negligible. For the case of the CP-odd Higgs boson also the full one-loop corrections 
to the decay into sfermions were analyzed in [5J E] ■ 

This paper is the continuation of the effort in [5J, [6J and includes the decays of the remain- 
ing Higgs bosons of the MSSM into sfermions (including crossed channels f2 — > fih°). 
It also extends the SUSY-QCD one-loop analysis of [7] by including all SUSY-QCD and 
electroweak effects. The emphasis is put on the decay into 3rd generation sfermions as 
their masses can be light due to large mixings. Nevertheless, analytical and numerical 
results are presented for all generations of sfermions (i.e. hi — > fifj and if ± — > where 
hi = (h°, H°) and / = (u, d, s, c, b, t, e, ft, f). 

The full electroweak corrections are calculated in the on-shell scheme [8] in the MSSM 
with real parameters. Due to the known problems of the on-shell scheme as demonstrated 
in [B], the artificially large on-shell parameters are replaced by the corresponding DR 
counterparts. The numerical analysis is made using the DR input defined by the Super- 
symmetric Parameter Analysis Project (SPA) [9]. In contrast to [6J, the actual calculation 
uses an on-shell input set fully consistent with the SPA convention. In order to obtain 
such a input set, the renormalization of the whole MSSM is required. 
The paper is organized as follows. In section 2 the tree-level formulae are given for all de- 
cays. Section 3 and 4 show the full electroweak corrections including the bremsstrahlung 
using the analytic formulae from the appendices A and B. The numerical analysis is 
presented in section 5 and section 6 summarizes our conclusions. 
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2 Tree-level result 

The tree-level widths for a neutral Higgs h® 12 y = {h°,H } decaying into two scalar 
fermions, h? k — > fifj with i,j = (1,2), are given by 

NL K,(mlo , m 2 ? , mf ) 

ptree/r _^ ff \ _ _ v n k h Si' \ r if |2 /i\ 

r V>, ■ ./../}) - j-^r I^W I 1 ) 

with k(x, y, z) = yj(x — y — z) 2 — 4y,z and the colour factor iV^ = 3 for squarks and 
Nq = 1 for sleptons, respectively. 

Analogously, the decay width for the charged Higgs boson H + is given by 

^K(ra|, + ) mi,ml) 

r tree (# + -////) = 77-^ — -^\g\i\\ (2) 

J lb7rm^ + J 

where y'^ stand for the up-type or down-type sfermions. The sfermion- Higgs boson 
couplings G{j k and Gjjx, defined by the interaction lagrangian £ int = G{j k h k f*fj + 
Gj^ H + fj* as well as all couplings needed in this paper, are given in [6]. 

The sfermion mass matrix is diagonalized by a real 2x2 rotation matrix R ia with 
rotation angle 9j pUl E] , 



( m h 


1* H 


V rn 2 RL 


™rr J \ 



mf a f m f \ + / mf 



which relates the mass eigenstates fi, i — 1,2, (m; < m^) to the gauge eigenstates / a , 

a = L, R, by Jj = R{ a f a - The left- and right-handed and the left-right mixing entries in 
the mass matrix are given by 

m| £ = M 2 q L} + (If L -e f sm 2 9 w ) cos2/3mf + mj , (4) 
mj^ = M 2 ^ q + e/ sin 2 6>vy cos 2/3 mf + m 2 , (5) 

a/ = Ay - ^(tan/3)" 2 ^ . (6) 

Mq, M£, Mj)-, Mg and Mg are soft SUSY breaking masses, Aj is the trilinear scalar 
coupling parameter, \i the higgsino mass parameter, tan f3 = ^ is the ratio of the vacuum 
expectation values of the two neutral Higgs doublet states pUl ELL], If L denotes the third 
component of the weak isospin of the left-handed fermion /, ef the electric charge in 
terms of the elementary charge eo, and 9w is the Weinberg angle. 
The mass eigenvalues and the mixing angle in terms of primary parameters are 
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m 



mf + m f - =F./(m 2 - — m 2 - ) 2 + 4a 2 f m 2 f ) , (7) 
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cos9f = / y (0<# f ~<7r) 
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and the trilinear soft breaking parameter Af can be written as 



— (m 2 h - m~J sin 28 j + (i (tan (3)~ 2I f L . (9) 



A f = 

The mass of the sneutrino v T is given by m 2 T = M| + | m| cos 2/3. 

For the crossed channels, / 2 — > A^fc and /J ~~ > f\H + , the decay widths are 



r tree (/ 2 -»■ Afyu = — *» J1 3 J2 l^kn (io) 

lb tx m% 

h 
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3 One-loop Corrections 

The one-loop corrected (renormalized) amplitudes G{ T k n and G||f cn can be expressed as 
G{fk n = G f ijk + AG f ijk = G{ jk + 6G{j% + 5G{^ + SG^Z , (12) 

where 6G{j k ,8G{j k and SG{j k and the corresponding terms for the couplings to the 
charged Higgs boson stand for the vertex corrections, the wave-function corrections and 
the coupling counter-term corrections due to the shifts from the bare to the on-shell val- 
ues, respectively 

Throughout the paper we use the SUSY invariant dimensional reduction (DR) as a reg- 
ularization scheme. For convenience we perform the calculation in the 't Hooft-Feynman 
gauge, f = 1. 

The vertex corrections SG^ and 5G\ji come from the diagrams listed in Figs. [TS] and 

[TBI The analytic formulae are given in Appendix [Bl The wave-function corrections SG(j^ 
can be written as 

dZi'i Gi> jk + SZffj G{ jlk + 5ZH G{ jt , (14) 



SG{^- 



with the implicit summations over = 1,2. The wave-function renormalization 

constants are determined by imposing the on-shell renormalization conditions [8] 

SZl = - Re tll(m% ), i = 1, 2 , 



fi 
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= 2 2 2 Reng(w&), M=(U), fc^i. 
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(15) 



(16) 



The explicit forms of the off-diagonal Higgs boson and sfermion self-energies and their 
derivatives, H kl , 11^. and IT^, tit are given in Appendix [A] and in [5]. 
The coupling counter-term corrections which come from the shifting of the parameters in 
the lagrangian can be expressed as 



8G$ 



6Rf ■ gL. ■ (Rff + Rf ■ 8G? ■ (R f ') T + R f ■ G? ■ (Stf) 



T 



.(17) 



The counter term for the sfermion mixing angle, 89j, is determined such as to cancel the 
anti-hermitian part of the sfermion wave- function corrections [T21 [T5] . Analogously we fix 
the Higgs boson mixing angle a by means of 



8a = \ (8Z« - 8Zg) = 1 Re (uf 2 (m 2 H0 ) + 11* (m^ 



Using the relations 



6G(n -Gfc, ^ = 4x> (19) 



da J da 

and absorbing the counter terms for the mixing angles of the outer particles, 8a and 89j, 
into 8G{j™ yields the symmetric wave-function corrections 



<r Symm ' ) = \^Zl + 8ZQ G{ 3k + \ (Szf,, + 5Zf, j )G{ j , k + \ (8Xj> + 8Xj>) G £, . 

(20) 

Note that in this symmetrized form momentum-independent contributions from four- 
scalar couplings and tadpole shifts cancel out. 
The sum of wave-function and counter-term corrections then reads 



£Qf{w+c) gQf(vi,symm.) 



R'-SG^-iR'f] (21) 
The explicit forms of the counter terms 8G^ LR k for k = 1 , 2 are given by 



($ g lr,i)ii = -\^h f m f c a (^^ + ^^^j-g z rn z e f 8s 2 w s 

+gz mz (lf-e,sl,)s^ m + *5E + JL) , (22) 
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for the sfermion couplings to the Higgs boson h° and 

($ G LR2)u = -V2h f m f s a (^ + ^^)+g z m z e f 5swC a+ p 

\ rif m,f J 

-g z m z (I 3 f L -e f s 2 w )c a+3 (— + - t a+g 5p) , (25) 
1 V g z m z J 

(^L, 2 )i 2 = 6 -^(G{^) 12 -^=(5A f s a -5^c a ), (26) 
(^L, 2 ) 2 2 = -V2h f m f sJ 



h f m f 

2 f s 9z . 5s^ \ 

g z m z e f s w c a+f3 1 h -5 £ a +/3 5/3 (27) 
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for the couplings to H°. 

Analogously to the decays of the CP-even Higgs bosons, the sum of the wave-function 
and counter-term corrections of the charged Higgs boson can be expressed as 



5G]$ w+c) = 5G]ji W ' symm '^ + R p ■ 6G 1 l R1 ■ (R 11 ) 1 .. + 5G\j[ w ' nvv+t1 ^' (28) 

L J iy 

with the symmetrized wave-function corrections 

^, symm .) = 1 {sz p + Mfctf^ . l( ( ;z;;, . 5Z^)G]^ + Uz^G\\ x . (29) 
The single elements of the matrix corresponding to the variation with respect to the 

T I 

couplings, oGl R1 , are given explicitly as follows: 



?/ T . XnU . (Ttf^T] J- xn M{w,HW+HG) 



{5G LR1 ) n - n, f m f sp\—— + — — + — 1 + ri f m f cp 

\ fit. writ, S3 j 



5h h ^Vt. &c B 



gm w . nQ (5g 5m w 5tan/3 

pr- sin 2p[ 1 h cos 2p — 

V2 V <7 rriw tan p 



(30) 



(^l J Ril ) 12 = -7-^(^,1)12 + ^/ 1(^/^/3 + ^/1^/3 + 5/i C/3 + fiScp) (31) 
"A 

(^,1)21 = -7-^(^.1)21 + h h {5A f ,cp + A fr 8cp + 5[iSf3 + fi5s B ) (32) 

(5C U ) h m r ( 5kf 1 4 4 ^ 4 h m « 4- ^ 4- 

(^ w ) 22 - / VT™/^( a ^ + ^ + — J+VW^"- + — + — J(33) 

The counter terms appearing in eqs. (I22H33I) can be fixed in the following manner. Some 
of them can be decomposed further as is the case for 5hf and 5g 

5hf 5g 5m f 5m w f — cos 2 (5 1 5 tan (3 5g 5e 5sin^ 

hf g mf mw \ sin 2 /3 J tan/3 ' g e sin6w 
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for | d "^ n }-type sfermions. 

For the remaining counter terms we use the standard renormalization conditions. The fix- 
ing of the angle (3 is performed using the condition that the renormalized A°-Z° transition 

2 

A" 



vanishes at p 2 = m 2 A0 as in |14j . which gives the counter term 



^ tan f = ^— M A o z »Ko). (35) 

tan [3 mz sin 2p 

The higgsino mass parameter \x is fixed in the chargino sector by the chargino mass matrix, 
<5/i = SX 2 2, as explained in detail in [15], [16]. 

The counter term to the Standard Model parameter sin^ is determined using the on- 
shell masses of the gauge bosons as in [T7|. To avoid the problems with light quarks in 
the fine structure constant a, we use the MS value at the Z-pole with the counter term 
given in [5], [21] . 

The on-shell counter term that has the biggest influence and also poses a serious problem 
is the counter term to the trilinear scalar coupling parameter Af. The explicit form of 
the counter term was already given in [5] and it was shown in E] that this counter term 
becomes very large for large values of tan p. One of the aims of this paper is to show that 
this problem is present in all Higgs decays into sfermions. The solution takes advantage of 
the fact that the SPA convention which we use here, defines the SUSY parameters in the 
DR scheme. Therefore, the trilinear scalar coupling parameters Af are taken DR without 
the use of the large on-shell counter term. 



4 Infrared divergences 

To cancel infrared divergences we introduce a small photon mass I and include the real 
photon emission processes h® — > fi fj'j {h\ = {h°,H }) and H + — > fjfj'y- The decay 
width of H+(p) fj(h) + fj(-k 2 ) + 7 (fc 3 ) (Fig. m> is given by 

it it 
(a) / h (b) / Ji (c) / h 



/ 

/ 



XA/^ 7 / 
H +S^L^ / H + < H + _, x 

STL ^ N ^ 

Figure 1 : Real Bremsstrahlung diagrams relevant to cancel the IR-divergences in H + (p) - 
fl{ki) + fH—k-z) + 7(^3). The diagrams for the neutral Higgs decays are analogous. 



T(H + - fj # 7 ) = N 3 C |Ggil 2 (-eo) 2 [m^ + Joo + e 2 m 2 / u + e 2 m 2 / 22 
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-e t e b ((m 2 H+ - m 2 - m 2 )I 12 - h - h) - e t ((m 2 - m 2 H+ - m 2 )I 01 - h - Io) 



+ e b ( (mi - m H+ - - h~ h 



with the phase-space integrals I n and I mn defined as [18 
cPki d 3 k 2 d 3 k s 4 



Ih - in ~ J m m m 5 (p ~ kl ~ h ~ h) (2^ + >F) . . . (2^ + ^ ■ (36) 

The full IR-finite one-loop corrected decay width for the physical value I = is then given 
by 

T co "(H + - f} ft) = r(H + ^fJf] l ) + T(H + ^fJlh) (37) 

Analogously, for the neutral Higgs boson decays and the crossed channels the photon 
emission processes yield 



('- 



(e oe f ) 2 \G{ jk \ 2 
167r 3 m^o 

2 



m h o -rrif -m f )I 12 + mJn + m? hi—I\—h 



nh^hKi) = ^^J^f [K- m l ~ m / 2 ) /o1 - m l In - m l Io ° - Io - h 

(38) 

where the IR-finite decay widths are 

T c ™(h° k -> Ul) = T(h° k -> + T(h° k -> 7 ) , (39) 
r corr (/ 2 - fih° k ) = T(f 2 -> M°) + T(/ 2 - M° 7 ) . (40) 



5 Numerical analysis 

The numerical results presented in this section are based on the SPSla' benchmark point 
as proposed by the Supersymmetric Parameter Analysis Project (SPA) [9]. A consistent 
implementation of the SPA convention into the calculation of a decay width and the nu- 
merical analysis is a non-trivial endeavor. As the electroweak one-loop calculations are 
carried out in the on-shell scheme and the SPA project proposes the SUSY input set in 
the DR scheme at the scale of 1 TeV, a conversion of the input values is necessary. This 
conversion requires the renormalization of the whole MSSM in order to transform the 
input parameters correctly. Moreover, the numerical analysis of a decay makes varying 
fundamental SUSY parameters necessary. That is why the above mentioned transforma- 
tion of parameters has to be carried out for every single parameter point. In our case 
this is provided by the not-yet-public routine DRbar20S which couples to the spectrum 
calculator SPheno [19]. The transformation is performed in the following two steps: 



S 



1. The SPA input, i.e. the on-shell electroweak SM parameters, the strong coupling 
constant and the masses of the light quarks defined in the MS scheme and the masses 
of the leptons and the top quark defined as pole masses and the SUSY parameters 
defined in the DR scheme at 1 TeV, is given to SPheno which transforms it to a 
pure DR input set including also higher loop corrections. 

2. The pure DR set is taken as input for the DRbar20S routine which yields as output 
the complete set in the on-shell scheme. An example of different sets of parameters 
for the SPSla' benchmark point can be seen in Table [TJ 

All plots below show the dependence of the decay width on a DR parameter. By varying 
a single DR parameter and transforming subsequently to the on-shell scheme, almost all 
parameters are influenced. That means, not only the corresponding on-shell parameter 
changes, but also the other parameters through loop effects. 



DR SUSY Parameters 


g' 


0.36354 


Mi 


103.21 


9 


0.64804 


M 2 


193.29 


9s 


1.08412 


M 3 


572.33 


Y T 


0.10349 


A T 


-445.4 


Y t 


0.89840 


A t 


-532.3 


Y b 


0.13548 


A b 


-938.9 




401.63 


tan (5 


10.0 


M Ll 


1.8121 • 10 2 


M L , 


1.7945 • 10 2 


M El 


1.1572 • 10 2 


M E , 


1.1002 • 10 2 


M Ql 


5.2628 • 10 2 


Mq, 


4.7091 • 10 2 


M Vl 


5.0767 • 10 2 


Mu, 


3.8532 • 10 2 


M Dl 


5.0549 • 10 2 


M D , 


5.0137 • 10 2 




2.5605 • 10 4 




-14.725 • 10 4 



On-shell SUSY Parameters 


9' 


0.35565 


Mi 


100.31 


9 


0.66547 


M 2 


197.01 


9s 


1.08419 


M 3 


612.81 


Y T 


0.10771 


A T 


-394.2 


Y t 


1.04638 


At 


-495.0 


Y b 


0.20481 


A b 


1197.8 




398.71 


tan (5 


10.31 


M Ll 


1.8394 • 10 2 


Ml, 


1.8199 • 10 2 


M El 


1.1784- 10 2 


Me, 


1.1172 • 10 2 


M Ql 


5.6390 • 10 2 


Mq, 


5.0369 • 10 2 


M Vl 


5.4540 • 10 2 


Mu, 


4.1021 • 10 2 


M Dl 


5.4352 • 10 2 


Md, 


5.3894 • 10 2 


M 2 Hl 


2.7220 • 10 4 


Mk 


-15.726 • 10 4 



Table 1: The input parameters for the SPSla' point according to the SPA project. 
Left: DR input values at Q — 1 TeV, Right: on-shell values used in the calculation. 

Comparing the parameter sets in Table [1] one can easily see that the on-shell counter 
term for A b is large as there is a huge difference between the DR and the on-shell value 
of the trilinear scalar coupling parameter. This is caused by fixing the of A b parameter 
in the sfermion sector [5]. This fixing was shown to lead to a numerically large counter 
term which should be avoided. The decays of Higgs bosons into sfermions (and the cor- 
responding crossed channels) are the only 2-body decays that are affected directly as the 
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trilinear coupling parameter appears at tree-level. In this case, a very large counter term 
makes the perturbative expansion unreliable. Here we make use of the fact that the in- 
put parameters are given in the DR scheme. That means our calculation uses on-shell 
parameters except for the Ab and m& which take the original DR values. Although not 
shown in the Table [U this behaviour is common to all down-type trilinear scalar coupling 
parameters, and the same strategy as described for the Ab is applied for them as well. 
A distinct feature of all decay modes involving down-type sfermions is a large difference 
between the on-shell and the SPA tree-level. The origin of this difference is again the 
large counter term for the trilinear scalar couplings. 

Keeping the numerical analysis strictly confined to the SPSla' benchmark scenario would 
mean that most of the possible decays are kinematically not allowed. The vertical red line 
denotes the position of the SPSla' parameter point for the kinematically allowed decay 
modes. For the other decays we slightly deviate from the SPSla' point adjusting mainly 
m^o and the relevant soft supersymmetry breaking terms Mjq fj £> gm- These parameters 
only influence the kinematics and have usually no effect on the couplings. 
In general, we always show the results using the on-shell parameters (dotted curve for on- 
shell tree-level and red dashed curve for on-shell full one loop decay width) as well as the 
improved decay widths where the parameters Af and rrif are taken DR (dash-dotted curve 
denotes SPA tree-level and blue solid curve stands for the full one loop decay width). This 
convention does not apply in cases where there is no down-type trilinear scalar coupling 
entering the tree-level. There we show only the on-shell and SPA tree-level together with 
the final one-loop decay width. For comparison with other calculations, the SPA tree-level 
is shown as defined in [9] taking all parameters in the couplings in the DR scheme and 
using the proper masses for the kinematics. 



{M D = 150 GeV, m 2 A0 = 10 6 GeV} 



^ 0.15 

o 

"3 0.05 
T -0.05 

o 

U -0.15 

5 10 15 20 25 

tan (5 

Figure 2: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 

tree-level (dash-dotted line) and full one- loop corrected width (solid line) of H° — > b±bi 
as a function of tan (3 according to the SPA convention. 
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Figs. El [3] and H] show the dependence of the decay width on tan /3 and clearly demonstrate 
the known fact [7J [5] that the counter term for Af, grows with tan (3. As mentioned above 
such a large counter term causes the perturbation series to break down. The full one-loop 
results in the on-shell scheme (red dashed curve) differ from those where DR parameters 
are used (blue solid) only by higher orders. Nevertheless, due to the perturbation series 
breakdown, the higher order corrections are no longer suppressed by the coupling con- 
stant. This is the reason why the results using the DR parameters should be viewed as 
the final one-loop corrected decay widths for the processes calculated in this paper. 
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Figure 3: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one-loop corrected width (solid line) of H° — > d\di 
as a function of tan f3 according to the SPA convention. 
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Figure 4: On-shell tree- level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one-loop corrected width (solid line) of H° — > fx fx 
as a function of tan (5 according to the SPA convention. 
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The next class of plots shows the dependence on the superpotential parameter /j, (see 
Figs. [3, El and [5]). The typical behaviour of the tree- level is governed by the square 

of the coupling G{j k . It implies that the tree- level is a quadratic function of n and as 
all one-loop corrections in this case are factorizable this dependence is preserved in the 
full one-loop result. In case down- type sfermions are involved, the on-shell curves are 
deformed by the large difference in the Af parameter. The corrections can reach up to 
40% for some areas of parameter space and are therefore not negligible. 

{M b = 150 GeV, m 2 A0 = 10 6 GeV} 



[GeV] 


3 - 


2 - 




1 - 




T 
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- 




-1 - 




Figure 5: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one-loop corrected width (solid line) of H + —* t\bi 
as a function of tan f3 according to the SPA convention. 



{M D = 150 GeV, m 2 A0 = 10 6 GeV} 




H [GeV] 

Figure 6: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one- loop corrected width (solid line) of H° — > S1S2 
function of n according to the SPA convention. 
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Furthermore, the pseudo-threshold in Fig. [7] comes from the sbottom contribution to the 
Higgs wave-function correction. 

The one-loop width of H + — > z> r f 2 (see Fig. [S]) is unexpectedly sensitive to the large 
difference of the on-shell and DR Ab parameters above fi = 600 GeV although there is 
no such parameter at tree-level. It is caused by the enhanced contribution of the vertex 
diagram with a 4-sfermion coupling. This diagram contains the coupling of the charged 
Higgs boson and a stop-sbottom pair where the Af, parameter appears. 

{Mfj = 150 GeV, m 2 A0 = 10 6 GeV} 




-1000 -500 500 1000 

fi [GeV] 

Figure 7: On-shell tree- level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one-loop corrected width (solid line) of H° — > tit 2 as 
a function of according to the SPA convention. 
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Figure 8: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one-loop corrected width (solid line) of H + — > v T i~2 
function of \i according to the SPA convention. 
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Fig. [9] illustrates the aforementioned problems of the perturbation series in case of using 
the on-shell Ab parameter. As one can see there is no obvious divergence of the decay 
width in the on-shell scheme. Nevertheless, the full one-loop widths in the on-shell scheme 
and the one with Af, taken DR are far apart. This separation is a pure two loop effect 
coming from using different Af, values when calculating the SAf, counter term. 
In this particular case the electroweak corrections interfere destructively with the QCD 
corrections reducing them by half. 
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Figure 9: On-shell tree-level (dotted line), full on-shell one-loop decay width (dashed line), 
tree-level (dash-dotted line) and full one- loop corrected width (solid line) of hi —* bih° as 
a function of ji according to the SPA convention. 




1 3 



£ °- 8 V 

a ■ > 

^ 0.6 

'■^ 0.4 - 

T 

0.2 - 



-1000 -500 500 1000 

A t [GeV] 

Figure 10: On-shell tree-level (dotted line), tree-level (dash-dotted line) and full one- 
loop corrected width (solid line) of t% — > t\hP as a function of A t according to the SPA 
convention. 
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The only plot over the trilinear coupling Aj is shown in Fig. [10] which is for the decay 
t 2 -> t x h Q . The region A t = (-120 GeV, 320 GeV) is kinematically forbidden. Although 
at some regions of parameter space (e.g. A t = (700 GeV, 900 GeV)) the correction to the 
SPA tree-level is large, one can see there is also a large difference between SPA tree-level 
and the tree-level used in the perturbation expansion denoted by the dotted line. 

6 Conclusions 

We have calculated the full electroweak one-loop corrections to the charged and CP-even 
neutral Higgs boson decays to sfermions including the crossed channels. We have also 
included the SUSY-QCD corrections which were calculated in [7j. Similar to [5] and [6], 
the on-shell parameters Aj, and m& (and the corresponding down-type parameters for the 
first two generations) were replaced by their DR values to avoid the numerically large 
counter term. Furthermore, we have presented the first consistent numerical analysis for 
a one-loop decay width based on the Supersymmetric Parameter Analysis project [9]. 
This required the renormalization of the whole MSSM in a way that allows to carry out a 
transformation between the on-shell and the DR scheme for every single parameter point. 
The corrections were shown not to be negligible and were comparable in the magnitude 
to the QCD in some regions of parameter space. 
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A Self-energies and counter terms 

Here we give the explicit form of the self-energies needed for the computation of the one- 
loop decay widths {h°, H } -> fj d and H+ -> f]fj. 

For the neutral and charged Higgs fields we use the notation h\ = {h°,H }, = 
{h°, H°, A°, G }, H£ = {H+, G + ,H-,G-} and = (H+)* = {H~, G~, H+, 
The sfermion self-energies and their derivatives as well as the vector boson self-energies 
and the scalar- vector mixing of A°Z° have already been published and can be found in 
|6J. Also most of the couplings used in this paper are given therein. 

In the following we use the standard two- and three-point functions Bi and Cj from [20] 
in the conventions of |18j . 



A.l Diagonal Wave-function corrections — derivatives of Higgs 
boson self-energies 

The conventional on-shell renormalization conditions for the diagonal wave-function renor- 
malization constants are given in terms of the derivatives of the corresponding self- 
energies, 



8Z% = -Renj£(m*o) 



(41) 



where the dot in IlJ^k 2 ) denotes the derivative with respect to k 2 . In the following we 
list the single contributions of the Higgs wave-function corrections. The derivatives of the 
CP-even Higgs bosons h° and H° depicted in Fig. [TT] are given as follows: 
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Figure 11: Diagonal self-energies of CP-even Higgs bosons h° and H c 
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where we have used 

Ml = 



■wx[a + /3- f(fc- 1)] cos[a + /3- §(&- 1)] 
cos[a + /3 - f (A - 1)] sinja + /? - §(jfc - 1)] 



17 



A'{k) A.(k) fa _^ _n\ 

mn mnKr- 1 HJ 
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\ sin 2/5 -cos 2/5 J ' 

The diagrams in Fig. [12] show the diagonal charged Higgs boson self-energies entering in 
the wave-function corrections. 
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Figure 12: Self-energy diagrams of the charged Higgs boson H + relevant for the calculation 
of diagonal wave-function corrections. 
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A. 2 OfF-diagonal Wave-function corrections — Mixing of CP- 
even Higgs bosons 

According to Section [31 the off-diagonal wave-function renormalization constants of the 
external Higgs bosons h° k = {h°, H } are given by 
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Figure 13: Diagrams for off-diagonal mixing of the CP-even Higgs bosons h° and H G 
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A. 3 # + VF + -mixing 



The scalar- vector mixing self-energy, Hjjw (k 2 ) (see Fig. [T4j) . is defined by the two-point 
function 
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B Vertex corrections 



B.l hlfifj vertex 

Here we give the explicit form of the electroweak contributions to the vertex corrections 
which are depicted in Fig. |T5l For SUSY-QCD contributions we refer to [Tj. 
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Figure 14: H + W + -mixing self-energies 



+ SG{$ VVS) + SG^ SS) + 5G{^ VV) + 8G$ hH ^ + 5G{^ (70) 

The single contributions correspond to the diagrams with three scalar particles (5G{j^,' H ^^ 
and 5G{p. ,HH ^), three fermions (5G{j%'* and SGf^'^^), three particles with one or two 



vector bosons (SG^ and 8Gjj k ' ) and two scalar or two vector particles {8G\j k 

r ijfc 



and 8G{^' VV ^) in the loop. The vertex corrections due to the mixing of the outer particles, 



i. e. Higgs and sfermion mixing terms 8G{j k ' hHm and 8G(j^'^ m will be combined with 
the counter terms of the Higgs and sfermion mixing angles, 5a and 89 z, see eq. ( )20l) . The 
vertex corrections from the exchange of one Higgs boson and two sfermions are 
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with the standard two-point function Co [2U] for which we follow the conventions of 
The graph with 2 Higgs particles and one sfermion in the loop leads to 
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For the gaugino exchange contributions we get 
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with the abbreviations h = (^oflW + 9o9ij2),hi = (go 9^2 + 9o9i92) and /i 2 = 
(9o9i92 + 9o9\9k)- For up- type sfermions F ; + fc = F,+ fc and for down- type sfermions 
chargino indices are interchanged, F^ fc = F+ ifc . 

We split the irreducible vertex graphs with one vector particle in the loop into the single 
contributions of the photon, the Z-boson and the VF-boson, 
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In order to regularize the infrared divergences we introduce a photon mass I. Thus we 
have 
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where we have used the vector vertex function 
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The vertex corrections coming from loops with two vector bosons and one sfermion are 
given by 
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We split the irreducible vertex graphs with two scalar particles into the contributions 
from two Higgs bosons, two sfermions and the corrections stemming from Higgs-sfermion 
loops. For better reading we introduce the abbreviations 
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All other combinations occurring are obvious. 
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The contributions due to the exchange of sfermions from the other two generations, F n 
are given by 
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where F denotes values belonging to the scalar fermions with the same isospin as /, but 
from the other two generations (e.g. F\ = {u\, c\} for the stop case, . . . ), and F' sfermions 
with different isospin etc. 
With the abbreviations 



'hi 



- kl 



/ cos a cos (3 cos a sin (3 \ 

sin a cos (3 sin a sin (3 

-zcos 2 /9 -f sin 2(3 

V -f sin 2/3 -isin 2 (3 J 



I — sin a sin /3 sin a cos /3 \ 

cos a sin (3 — cos a cos (3 

-isin 2 /3 | sin 2(3 

V i sin 2(3 -i cos 2 /? / 



26 



Jb,0+ 
-kl 



/ — sin a cos (3 — sin a sin (3 





\ 


cos a cos (3 cos a sin f3 










-f sin 


2(3 -i sin 2 (3 










\ z cos 2 


13 | sin 2/3 





o J 


/ 


cos a sin /3 — cos a cos 









sin a sin /? — sin a cos 


p 







— | sin 2/3 zcos 


2 (3 






V 


— i sin 2 /3 | sin 


2(3 


J 





Jt,0+ 
-kl 



the diagrams with one Higgs boson and one sfermion in the loop lead to 
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with /if = {/i/ T , 0} and = hf } for the decay into {squarks, sleptons}, respectively, 
z «-> j stands for both terms before with % and j interchanged. 
Finally, for the vertex graphs with two vector bosons we obtain 
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B.2 tf+jj/j vertex 

In this chapter we list the various contributions to the vertex corrections of the charged 
Higgs boson decays H + —> ///i. For simplicity the formulae are given for third generation 
sfermions as the substitution for the first and second generation sfermions are obvious. 
Using the definitions for various generic vertex functions and products of couplings from 
the previous chapters we get for the vertex corrections, 

SG fb{v) = SG^ V,H ^ + SG fb{v,HHf) + 5G fb(v '* ff) + 5G tb{vjm + 8G fb{v ' lSS) 

ijl ijl ijl ijl ijk ijl 

+ 5G^ ZSS) + 5Gf^ wss) + 5G^ HH) + 5G^ m + 5G^ PP>) 

+ 5G^' Hf) (93) 
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The single contributions are given as follows: 
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Figure 15: Vertex diagrams relevant to the calculation of the virtual electroweak correc- 
tions to the decay width h° k -> fj j with h° k = {h°, H }. 
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Figure 16: Vertex diagrams relevant to the calculation of the virtual electroweak correc- 
tions to the decay width H + — > ///,-. In the fourth row, /J and denote up- and 
down-type sfermions of all three generations, respectively. 
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